arXiv:1508.05975vl [cond-mat.stat-mech] 24 Aug 2015 


Bound States of Majorana Fermions in 
Semi-classical Approximation 


G. Mussardo^’^’^ 

^SISSA and INFN, Sezione di Trieste, via Bonomea 265, T34136, Trieste, Italy 
^International Centre for Theoretical Physics (ICTP), 1-34151, Trieste, Italy 
^International Institute of Physics, Natal, Brasil 

We derive a semi-classical formula for computing the spectrum of bound states made of Majorana 
fermions in a generic non-integrable 2d quantum field theory with a set of degenerate vacua. We 
illustrate the application of the formula in a series of cases, including an asymmetric well potential 
where the spectra of bosons and fermions may have some curious features. We also discuss the 
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I. INTRODUCTION 

The experimental advances in neutron scattering and material design, as well as in the manipulation of 
ultra-cold atom systems, have recently permitted to realise and study in detail the class of universality of 
some of the most important low-dimensional quantum field theories. This is the case of the 2d Ising model 
in a magnetic field mm. a model not only important for its physical and historical relevance but also for 
its deep relation with beautiful mathematical object such as Es Lie Algebra. The class of universality of 
this model has been experimentally accomplished by designing a quantum spin chain made of CoNb 206 
(cobalt niobate) whose properties were subsequently analysed by neutron scattering [3] . The experimental 
analysis has confirmed, in particular, the theoretical prediction about the existence of three bound states 
below threshold [1], with their weights in the spectral density of spin-spin correlation function being in 
perfect in agreement with their theoretical determination [2]. Another example of class of universality is 
the celebrated Sine-Gordon model gHS], recently experimentally realised in terms of a cold-atom set-up 
and then studied through interference patterns [7]. This experimental work has permitted, in particular, 
to study the phase correlations, to characterise the topologically distinct vacua and the interactions of 
the excitations above them. 

Given the present trend of such a theoretical and experimental advance, one may hope to witness 
progress in the next future on one of the long-standing aspects of quantum field theories (QFT), alias 
how to determine the spectrum of their physical excitations. This is a genuine dynamical problem, since 
it is a well known that the QFT Lagrangians may employ fields whose particle excitations are actually 
not present in the physical spectrum. 

The most notorious case is probably Quantum Chromodynamics, where the quarks and gluons of its 
Lagrangian are absent from the spectrum which is given, instead, only by their hadronic and mesonic 
bound states (see, for instance, 0). Let’s also remind that for theories as QCD the only reliable way 
known up to now to get the masses of the hadron/meson particles is through numerical approaches - 
such as Montecarlo simulations of lattice gauge theories (see e.g. [3]). 

Another famous example - which has the advantage to be exactly solvable, therefore it neatly clarifies 
what kind of problem is to deal with the spectrum of a QFT - is the aforementioned Sine-Gordon model: 
it describes the dynamics of a two-dimensional bosonic field with interaction given by the potential 

Usg{t) = ^ [1 - cos(^¥>)] , (1) 

where mg is a mass-like parameter and /3 is a coupling constant. Such a theory has an infinite number 
of degenerate vacua | n) (n = 0, ±1, ±2,...), localised at tpn'^ = 27rn//3, each of them with the same 
curvature of the potential Usg{.t) these points, namely = m^. In terms of the usual and familiar 
perturbative arguments, one would be then tempted to conclude that each vacuum has always, at least, 
one neutral excitation above it. However this conclusion may be false. Indeed the exact ^-matrix of this 
model [4] shows that the situation may be rather different and even quite surprising: indeed, if /3^ > 47r 
such a particle - the very same particle described by the field ^p{x) itself - does not exist! When this 
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happens, as in the analogous case of QCD, one shall conclude that the spectrum of the theory is not made 
of the field/particle entering the Lagrangian but consists of other excitations that have to be properly 
identified. 

We brought these two examples simply to argue that the determination of the spectrum of an interactive 
QFT may be quite a difficult and non-perturbative problem which, apart from numerical methods, can be 
in general addressed by few analytic approaches as, for instance, the Bethe-Salpeter equation [10]. There 
are cases, however, where life is easier. This happens, for instance, for two-dimensional integrable QFT’s. 
In this case, indeed, the infinite number of conservation laws implies the elasticity and the factorization of 
their S'-matrix mm, which can be then exactly determined in innumerable cases (see, for instance, El); 
and, of course, when the exact expression of the S'-matrix is known, the spectrum of the corresponding 
theory tidely derives from a proper identification of its poles. 

Luckily enough, certain progress can be also made when the two-dimensional theories are not integrable, 
depending of course on the degree of ^^non-integrability", so to speak. For instance, if the theory in 
question may be regarded as small deformation of an integrable model, one can then apply the Form 
Factor Perturbation Theory (FFPT) [r^HTC] to follow the evolution of the spectrum moving away from 
the integrable direction. Eventually combining this method together with the efficient numerical method 
of Truncated Conformal Space Approach (TCSA) [T^, a big deal of information can be obtained on a 
series of aspects of QFT, such as confinement of topological excitations, adiabatic shift of the mass of 
non-topological excitations, decays of higher mass particles, presence of resonances etc. (see, for instance, 
[HHIlllIlHSSj). It is also worth mentioning that for gauge theories in two dimensions, the Bethe-Salpeter 
formalism finds an elegant application in the computation of the mass spectrum of the mesons |24j . 

In addition to all methods above mentioned (FFPT, TCSA and Bethe-Salpeter equation), there exists 
another one - the so-called semi-classical method |25j - which has the advantage to rely neither on 
integrability nor of small breaking thereof. While, in general, the application of this method may be 
quite intricated, it drastically simplifies if applied to the analysis of QFT’s with vacua degeneracy and 
topological kink excitations, as shown in |26H29j . In particular, properly refining and interpreting the 
relevant matrix elements on kink states originally stated in the paper by Goldstone and Jackiw |26j . 
in the paper [28] we have shown that it is possible to carry out a thorough study of the semi-classical 
spectrum of neutral bound states in purely bosonic kink-like theories. 

The aim of this paper is thus to extend the analysis of the paper |28j to the case in which there are 
also Majorana fermions, taking as our starting point the paper by Jackiw and Rebbi [29] . As we will 
see, there are a series of quite interesting situations that can analysed with such a formalism, such as: 
(a) QFT with the bosonic part of the potential having symmetric or asymmetrical potential wells; (b) 
QFT with kink excitations which are invariant under Supersymmetry; (c) QFT which are integrable or 
non-integrable. 

To make the analysis of this paper self-contained, in Section |ll| and |lll| we will briefly review the main 
results relative to the semi-classical treatment of bosonic theories. In Sect ion [TV] we address the topological 
excitations of QFT with fermions and the new vacua structure induced by them. Section]^ is the central 
part of the paper, where one can find our semi-classical formula for computing the spectrum of fermions 
in theories with degenerate ground states. The remaining sections are devoted to the analysis of a series 
of relevant examples, among which: (i) QFT with symmetric and asymmetric vacua (Sections VI and 
VII) (ii) Supersymmetric theories (Sections jVIII and IX); (hi) Integrable (Section]^ and non- integr able 


Supersymmetric theories (Sections jXIj and XII). Our conclusions are then discussed in Section XIII 


II. PURELY BOSONIC THEORIES 

Let’s initially consider a purely two-dimensional bosonic theory since it will set the stage for all our future 
considerations. Such a theory concerns with a scalar real field g}{x) and a Lagrangian density given by 

, ( 2 ) 

where A is some coupling constant. The case of our interest is when the potential I7((/?, {A}) possesses 
several degenerate minima, localises at (a = 1) 2,..., n), with U = 0 as shown in Figure 1. In 
QFT language, these minima identify different vacua | a) and therefore the basic excitations of such a 
theory are the topological configurations that interpolate between any neighbouring pair of them, namely 
kinks and anti-kinks. Let’s see in more details what they consist of. Semi-classically such excitations are 
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FIG. 1: Potential U{ifi) of a bosonic quantum field theory with kink excitations. 


identified as the static solutions of the equation of motion, alias 


with boundary conditions 


dl^ix) = C/'[^,{A}] , 


V3(-Oo) = 
ifi+OO) = 


& = a ± 1 . 


(3) 

(4) 


Let’s (fiabix) be the classical solutions which interpolate between (at x = —oo) and (at x = +oo): 
with 6 = 0+1, Pabix) is conventionally called kink while Pba{x) anti-kink. In the following, however, we 
will often find easier to use the word “kink ”to denote both of them, especially in discussions of general 
context. It is worth stressing that these kink configurations can be equivalently obtained as solutions of 
the simpler first order differential equation (with the same boundary conditions as above) 

^ = ±^/^ , (5) 

where the ± signs refer to the kink and the anti-kink respectively. 


Basic properties of the kink/anti-kink solutions. It is useful to remind some features of the 
kink/anti-kink solutions particularly relevant for the physical consequences they induce. In the following, 
unless differently stated, 6 = a + I. 


I. As it will become soon evident, it is important to know how pab(x) approaches the two asymptotic 
values p^^'^ (i = a, 6). This problem was solved in [28], where it was shown that when the curvature 

uji = U''{pf^) {i = a, b) at the two vacua is different from zero, these approaches happen in an ex¬ 

ponential way. The simplest way to show this result is to introduce the variables Tji = ^p(x) — p^^^^ 
{i = a, b) and expand the right-hand side of eq. ([^ (for our choice of the indices a and 6, the one 
with positive sign) in power of rji in the vicinity of the minimum p^'^ (for x —)■ —oo) or (for 
X —)■ +oo) so that 


drji 

dx 






= + a^'7?, -\-ayr]^ + . 


( 6 ) 


The solution of these differential equations provides the sought asymptotic approaches to the two 


vacua 


¥>(cr) = + , cr ^-oo (7) 

n—2 
oo 

p{x) = + 

n—2 


X —>■ +00 


( 8 ) 
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FIG. 2: Kink and antikink configurations. 


where the coefficients ^n^’s can be iteratively computed in terms of the an'^’s in Notice that 
all the exponents of the exponentials are expressed in terms of integer multiples of the curvature 
uii of the corresponding vacuum. When to a ^ toi,, the exponential approaches to the two vacua are 
of course different and in this case the kink configuration Lpab{x) may be not related in a simple 
way to the anti-kink configuration ipbaix)-, say by a parity transformation. As discussed later in this 
paper, this may have far-reaching consequences on the spectrum of the theory. 

2. Associated to each solution ipab{x) there is a classical energy density given by 

(^ab{x) = ^ +U{g}ab{x)) , (9) 

which, in light of ([^, can be written as 

eab{x) = 2U{(pab{x)) ■ (10) 


Integrating this density on the entire real axis, we get the classical value of the masses of this 
excitation 


M„h = 


Cabix) dx 


( 11 ) 


Notice that one does not need to know the exact solution of the kinks to get these values. Indeed, the 
masses Mat can be computed in terms of the potential U(ip) alone: using the monotonic behavior 
of the kink solution, making a change of variable t = fabix) and using once again eq. (HI , we have 
in fact 

Vb 

Mab = J v/2 u{p>) dip . (12) 

V>a 

Generally the masses Mab of the kinks is a non-perturbative quantity of the coupling constant A of 
the theory, which diverges when A —>■ 0. Hereafter a semi-classical regime has meant to be as the 
regime where the coupling constant is sufficiently small, so that the masses of the kinks are higher 
than any other mass scale. Eq. (El shows that, even though the kink (pab{x) and the anti-kink 
ipbaix) may not be related one to the other in a simple way, their mass is nevertheless always the 
same. 

3. Using the relativistic invariance of lagrangian theory (H. we can set the static solutions (pab{x) 
in motion by a Lorentz transformation ipab(x) —>■ tpab [(a^ ± vt) j^/l — : these new configurations 

correspond to lumps of energy fulfilling a relativistic dispersion relation and localised at the maxi¬ 
mum of the energy density eab{x). We can then associate to them the quantum kink states \ Kab{0)) 
HSHmiS], where 9 is the rapidity variable which parameterises their relativistic dispersion relation 

(13) 


E = Mab cosh0 


P = Mab sinh0 . 
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FIG. 3: Multi-kink configurations. 


At the quantum level, the masses Mab of the various kinks will experience a finite renormalization 
with respect to their classical values (121 (see [25l[3T]) and we will take explicitly into account this 
fact in some of the examples discussed in the sequel of the paper. 


4. The various kink configurations can be conveniently represented in a graphical way as in Figure 
Even though this graphical representation is a very crude image of the actual exponential approaches 
to the two vacua (which, as we saw, could be also different), it provides nevertheless a useful rule 
of thumb to build up and depict the multi-kink states of the quantum theory: these are given by 
a string of kinks, satisfying the adjacency condition of the consecutive indices for the continuity of 
the field configuration, as shown in Figure 

\ Ka,,aM Ka,,a,{e2) Ka,,aA93) ■ ■ ■) , (a,+i = tt, ± 1) (14) 

Notice that the field theory ([^ has an associated conserved current j^(x) = M whose 
charge 


OO 

Q = Af J di(p{x) = A/’[(/7(-I-oo) - (/7(-oo))] , 

— OO 


(15) 


assumes integer number values, provided we have properly chosen the constant Af in reference to 
the various values of the vacua ipn'^. These integers Q label the different topological sectors of the 
theory, where Q = Q refers to the Hilbert space build upon the various vacua | a), while <5 = ^1 to 
the sectors of kink and anti-kink. Other values of Q refer to higher multi-kink sectors [3T] . 


5. Important dynamical data are the bound states made of the kink-antikink excitations, alias the 
scalar excitations | Sc(^))aj called breathers, which may exist on top of each vacua | a). These 
excitations live in the Q = 0 sector of the theory. For a QFT based on a single real field, these 
breather are all non-degenerate states and, in presence of a charge conjugation symmetry C, they 
may be also classified in terms of the even or odd eigenvalues of C. On a general ground it may 
be argued that all scalar excitations that live on top of the vacuum | a) must be identified as the 
bound states of the kink-antikink configurations that start and end at the same vacuum | a), i.e. 

I Kab{9i) Kba{92) )■ These bound states correspond to the “tooth” configurations shown in Figure]^ 
Notice that if the vacuum | a) is connected to two different vacua rather than only one, there are two 
different kink-antikink states starting and ending at | a), alias | Ka^a±iKa±i,a)- However, as clarified 
in [3H], the only one that matters for the spectrum of the breathers is the configuration with the 
lower kink mass between the two values Ma,a±i or, in presence of a degeneracy of | Ka^a±iKa±i,a), 
a special linar combination of them. This is the content of the rule of “T/ie importance of being 
smair that we will also advocate later in this paper. So, let’s assume we have selected the kink- 
antikink configuration | Kab{ 9 i)Kba{d 2 )) associated to the lowest mass of the kink going out 
or coming to the vacuum | a): we can then identify the breathers | Be )a in terms of the poles at an 
imaginary value i ^ within the physical strip 0 < Im 0 < tt of the rapidity difference 9 = 9i — 62 
of this state 

9lb 


I Kab[9i) Kba{92)) 


I Be)a 


(16) 
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where gJJj, is the on-shell 3-particle coupling between the kinks and the neutral particle. Knowing 
the resonance value iu'^h, the mass of the bound states are simply obtained by substituing this 
quantity into the expression of the Mandelstam variable s of the two-kink channel 

s = 4:M*^ cosh^ ^ ^ TOc = ±2M*f, cos ^ . (17) 

Obviously for a bosonic spectrum we have to choose above the -I- sign: this will not be necessarily 
the case for the fermionic case, as discussed later. The determination of these resonance values 
u'^Ij and the semi-classical computation of the spectrum of these bosonic neutral excitations are 
discussed in the next section. 


III. THE BOSONIC BOUND STATES IN Q = 0 SECTOR 


When the coupling constant A goes to zero, the mass of the various kinks becomes very large: in this 
case there is a semi-classical self-consistent way to compute the matrix elements on kink states of the 
fundamental field f{x). The final formula just employs the Fourier transform of the classical field ^Pabix) 
interpolating between the two vacua | a) and | b) and reads 

OO 

= {Kabie2) I <7’(0) I KabiOi)) ~ J ipabix) , (18) 


where 0 = 9i — 02- This formula can be also generalised to any operator G[ip{x)], function of ipi^x): 
indeed, one has simply to substitute above (pab(x) G[ipab{x)]. The justification and the derivation of 
this formula are given in the original papers [2SH55] (see, in particular, [55] for its correct dependence on 
the relativistic invariants of the two-body channel) and will not further discussed here. 


One can now take advantage of such a formula to determine the semi-classical spectrum of the bound 
states in the Q = 0. As we argued in |55|, the steps to be taken are the following. Firstly, go to the 
crossed channel of the matrix element by substituting in eq.(18) 0 ^ in — 0, so that the corresponding 
expression may be interpreted as the Form Factor 


KbiO) = = {a I <p(0) I Kab{0i)KM)) , (19) 

where it appears the <5 = 0 neutral kink states around the vacuum | a) that we are interested in. Secondly, 
identify the poles of F^i,i0) localised in the physical strip 0 < Im0 < tt, see Figure]^ Since the poles of a 
Fourier transform of a function are given in terms of the exponential asymptotic bdiavior of the function 
itself, to get the spectrum above the vacuum | a) we have then simply to know the exponential approach 
of (pabix) to (pi^^ for X — >■ — OO. This asymptotic behavior is in eq. given that it consists of a series 
made of multiples of the same exponential (expressed by the curvature Wq at that minimum), the poles 
of F^ij{0) are regularly spaced by 




Ujg 

'xMab 


( 20 ) 


Focusing the attention on the kink of lower mass starting from the vacuum | a), the semi-classical mass 
spectrum of the bosonic neutral bound states around the vacuum | a) assumes then the universal form 


mf(n) = 2M:, 



= l,2,...Af 


( 21 ) 


la+l> 
I a> 

I a-l> 


FIG. 4: Kink-antikink configurations which may give rise to a bound state nearby the vacuum | a). 
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where is the total number of bound state on top of it: this is the number of poles that fall within the 
physical strip and it is given by 


= f- 

. sa. 

where [a:] expresses the integer part of the number x. Posing 


(1) = 2Mab sin ( 




the formula (211 can be eqnivalently written as 

sin (n 


(n) = 


sin ( 


n = l,2,...N^ 


( 22 ) 


(23) 


(24) 


Notice that, if > 1 the pole(s) of are outside the physical strip and there is no single bound state; 

in the two intervals ^ < 1 and | < , there are one and two bound states, respectively; while 

for < I, when the theory is non-integrable, out of the possible [l^a] bound states only the lowest two 
will be stable, the others being resonances [28]. Indeed, for a geometrical property of the sine of multiple 
angles, we have to„ > 2mi (for n > 3) and therefore these higher particles are generically expected to 
decay. Hence, the conclusion is that, in a generic non-integrable bosonic theory, each vacuum cannot have 
more than two stable neutral particles above it. 


IV. FERMION IN A BOSONIC BACKGROUND 

Let’s now consider a two-dimensional theory involving a scalar field ip{x) and a fermion field 

iKx) = ( * ) (25) 

with Lagrangian density 

^ -U{(p,X) -V{ip,g)4’'ip . (26) 

As coupling constants we have now both A and g. In addition to the bosonic potential I7((/?, A) (that we 
will take as those in Section 2, alias with a degenerate set of minima connected by the kinks ipabix)), we 
have now a Yukawa interaction V(tp) for the fermion. For this interaction we require that it satisfies the 
following conditions for any kink configuration ipab{x) 

Vi^M.g) = I 


X —>• —00 
X —>• +00 


(27) 




FIG. 5: Graphical representation of the pole and related bound state \ Be)a of kinks \ K^b) and anti-kink \ Kba) ■ 
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where both limits are finite and such that their product is always strictly negative 

i’a{g)vt{g) < 1 . (28) 

From now on, in order to simplify the notation, we will denote by Vab(x) the potential V{ip,g) computed 
on the kink configuration tpabix) and, for simplicity, we will also skip the presence of the coupling constant 
g in this expression. 

For the 7 matrices we adopt the Weyl representation given by 

7° = ^2=(®"o) ;7'=*^i=(.o) - 7^ = 7V=^3= (J . (29) 

The charge-conjugation matrix C, satisfying 

(C7°) (t'^)* (C'7°)-' = -r , 

is given by C = 7 *^ and maps the fermion ^ to its conjugate particle '0c according to 

0c = (C7°)0* . (30) 

In this representation a Majorana fermion (which satisfies the neutrality condition 0c = 0) bas then 
both components real 


01 = 01 


02 = 02 . 


(31) 


it was considered the situation where the Lagrangian (26) employs a complex Dirac 


While in 

fermion - and this leads to the discovery of the fractionalization of the fermion quantum number 
- here for simplicity we are interested in studying directly the case where ip{x) is a real Majorana 
fermion. The two situations are closely related, of course, and indeed our final result for the fermionic 


bound state spectrum in the <5 = 0 sector, eq. (58 1 , applies equally well to both cases. However, the key 


difference between Dirac/Majorana fermions consists of the impossibility to define, in the Majorana case, 
a fermion number as the conserved quantity associated to the transformation 0 —^ 0 , simply because 
the Majorana fermion is real. So, in the case of Majorana fermion, we will simply talk about a fermion 
parity V. A field theory similar to the one discussed in this paper was also considered in [3(7| . 

Notice that the Lagrangian (26) has no explicit fermion mass term and therefore the fermion looks 


massless. Yet the problem addressed in this paper is precisely how to find the fermionic mass spectrum 
in the Q = 0 topological sector. As a very preliminary observation in this direction, it is clear that it will 
be the interaction term V{(p{x)) to play the role of an effective mass for the fermion field 0(a:) and this 
is particularly true when the bosonic field ^p{x) has been prepared in a Q = 0 trivial state, such as one of 
its vacua: in these cases, the the effective mass of the fermion at the vacuum | a) is given by the constant 
value = Va- However, as we saw it was the case for the bosonic spectrum, to unveil the actual 

spectrum of the fermion in the Q — 0 sector (and see, in particular, whether the particle associated to 
tf{x) belongs or not to the physical spectrum!) we need to study first the higher topological sectors of 
the theory. 

Let’s consider then the Q = 1 sector (analogous considerations hold for Q = —1), with the bosonic 
field prepared in the one kink static configurations ipabix). Although, strictly speaking, such a classical 
configuration refers to the equation of motion @ where the fermion field is absent, yet it is natural to 
assume (as we will do) that the latter will have a negligible effect on the former as far as the mass of the 
kink itself far exceeds the “mass energy” of the fermion for both vacua, i.e. if 


Mab »| Uj I , i = a,b . (32) 

This will be the condition expressing the semi-classical limit of the Lagrangian ( |2^ . In the Q = I sector 
the fermion has a new qualitative non-perturbative feature: it possesses localised zero-energy mode, 
namely a normalizable solution of zero energy of the Dirac equation [25l[29]. Indeed, being ipab{x) a static 
configuration, we can as well as look for a static solution of the Dirac equations, which therefore take the 
form 


0a;02°^ + K 6 (a;) 0 f^ = 0 , 
9x0^^ + Vabix) 0^°^ = 0 . 


( 33 ) 
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FIG. 6: Kink configuration ipabix) and the associated fermion zero mode The latter is always localised at 

the point xg where the former makes a transition between the two vacua. 


These are nothing else but the 0-eigenvalue condition Hjj — 0 for the Dirac Hamiltonian 

Hd = -iadx + Vab{.x) 13 , (34) 


where a = 7 ^ 7 ^ = <73 and /3 = 7 *^. Solutions of the two equations (33) can be found in terms of the linear 
combinations of the two components of the fermion 




±0^°^) = A± exp 

X 

T J Vab{t)dt 


. Xq 


(35) 


In light of the asymptotic behaviors (28) of the potential Vab{t) and the condition (28) on their signs, it 
is clear that only one of the two solutions {x) is normalizable, the other being divergent at infinite. 
To fix the ideas, let’s assume that Ua < 0 and Vb > 0: in this case, the normalizable solution is 'ip+{x) and 
correspondingly we have to impose = 0 , so that the final result for wave function of the zero mode is 


'^abi^) = A tp exp 


X 

- J Vab{,t)dt 


(36) 


where '0 = ^ ^ ^ and H is a normalization constant. This wave function is usually localised at the 

position where the kink solution jumps between the two vacua, as shown in Figure 

In the Q — 1 sector, the quantization of the fermion field is achieved by its eigenmode expansion p51E5] 

0(a;, t) = ao (x) -)- ^ [bp 'tpE^ (x) -f bl ip-E^ (x)] , (37) 

P 

where the operator oq is relative to the zero-mode solution, while bp and are the operators associated 
to the positive/negative energy solutions '4>±e{x) of the Dirac Hamiltonian in the kink background (pab(x) 

Hd'P’±e{x) = ±Etp±E{x) . (38) 

The operators bp and bj^ satisfy the anti-commutation relations 

{bp,bl} = S{p-k) , {5p,ao}, = {6^,ao} = 0 . (39) 

The operator qq is instead a Majorana self-conjugated zero mode operator, oq = aj that, properly 
normalised, fulfills the condition 


( 40 ) 
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The ao mode carries zero energy and therefore the original quantum state | Kab{0)) has to be actually 
doubly degenerate [29]: we will use the notation | K^i^{0)) to identify these two degenerate states, related 
one to the other by 


I K+{9)) = ao I K-,{e)) , 
I = «o I K+{9)) . 


(41) 


Notice that we can consider | to be the charge-conjugated state of | and viceversa. Given 

their symmetry in sharing the Majorana fermion, there is a certain arbitrariness in assigning the fermion 
parity: our choice is regard | K^^^{9)) as kinks that have Majorana fermion parity V = ±1. The same 
considerations apply to Q = —1 sector, alias the sector of the anti-kinks | K^a), which now become 
doublets I K^^{9)). We will graphically distinguish these new kinks by putting on their previous graphical 
representation a red up arrow on the those of fermion parity -1-1 while a green down arrow on those of 
fermion number —1, as shown in Figure]^ 

We can now use this double multiplicity of kinks and antikinks to unfold the original vacuum structure. 
Namely, consider two bosonic neighborhood vacua | a) and | b) connected by the kinks and antikinks of 
different fermion number; then, the adjacency properties of the initial vacua, carried by the kinks | K^f), 
can be equivalently expressed as the adjacency conditions of a new set of vacua, in which we split one of 
the original vacua (say | a)) into two new ones | a±), while we leave the other | b) untouched, and we 
connect this new set of vacua as follows: | a+) is linked to | b) only by kink/antikink of fermion parity 
-1-1, while I a_) is connected to | b) only by those of fermion parity —1, see Figure]^ This construction 
can be repeated for any two neighborhood vacua and will be useful later. 


Let’s close this section by briefly discussing the non-zero energy solutions of the Dirac Hamiltonian 
equation (38) in the Q = 1 sector: notice that if 


i’E = 


( i’v 

V V'? 


stays for a positive energy solution, then the solution corresponding to the negative energy is simply 
given by 



a + 1) 
a) 

Kb) I K) 

a + 1) 
a) 

Ka) I Ka) 




FIG. 7: Graphical representation of the fermionic kinks \ and anti-kink \ The red up arrow denotes 

those of fermion parity V = -\-l, while the down green arrow those of fermion parity V = —1. 
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FIG. 8: (i) Adjacency diagrams of original vacua connected by fermionic kinks; (li) splitting the original vacuum 
I a) into two vacua \ a±), where each of them is related to the vacuum \ b) by kinks of the same Majorana fermion 
parity. 


alias there is always a pairing between the positive/energy Hilbert space. Moreover, the eigenvalue problem 
(381 can be equivalently expressed in terms of two Schrodinger equations 


-dl + (Ki + K'h) 


u = E'^u ; 

V = E'^ V , 


(42) 


for the linear combinations of the two components of these spinors 


U = Iprj+lp^ , 

V = ■ 

It is worth noticing that the equations ( |42[ ) are those of one-dimensional Supersymmetry Quantum 
Mechanics and shape-invariant potentials [331 134] . an observation that may be quit useful to find an 
explicit form of eigenvalues and eigenvectors in the Q = ±1 sectors (see, for instance, [3S]). 


V. THE FERMIONIC BOUND STATES IN Q = 0 SECTOR 

Let’s now study the fermionic bound states spectrum in the Q = 0 sector. To this aim, it is useful to 
remind that in the paper |29| Jackiw and Rebbi derived semi-classical matrix elements of the fermion 
field ifix) on the kink states | Actually, as they are written in [^, these matrix elements do not 

take into a proper account the dependence on the relativistic invariants of the channel of the two kinks. 
But this is a problem that can be easily cured, similarly to what was already done for the pure bosonic 
case in HZlllH]: one has simply to introduce the rapidity variable 9 and consider that, in semi-classical 
approximation, Mab sinh0 ~ Mabd, given that the kinks must be regarded to be almost static being their 
mass very large. In this way, the definite semi-classical formula for the fermion matrix elements on the 
kink states reads as 

OO 

5a.(^) = 9abi0) ^ J , (43) 

— OO 


where we have put 


9abi(^) = I m I K^b{02)) . 


(44) 
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In eq. (43l, 9 = {9i — 62 ) and is the zero-mode configuration of the fermion in the presence of the 


kink background Lpab{x). The identity of the two matrix elements is due to the real condition on 

the Majorana fermion ip{x). For a Dirac fermion field these matrix elements are complex conjugate one 
to the other and eq. (43) must be substituted instead by the identity 


{K^,{e,) I V'(O) I = {K^,{92) I V'^(O) I 


(45) 


Having now the correct dependence of this matrix element on the relativistic invariants of the channel of 
the two kinks, as in the bosonic we can make the analytic continuation 6 ^ m — 9 oi (43) and go to the 
crossed channel 


GUQ) = = {a I V'(O) I K^b{(^i)KU02)) ■ 


(46) 


Therefore, for the (5 = 0 bound states with fermion quantum number around the vacuum | a), we have 
to look at the poles of (?^^(0) in the physical strip 0 < Im0 < tt. Notice that there are two states which 
will share the same poles, alias | and | these states are charge-conjugated one to the 

other, in agreement with the general result that any solution with +E is paired with another of —E. 
Notice that, unfolding the original vacua, these states | have the graphical representation shown 

in Figure 


Given both the conditions (27) and (28), it is clear that zero-mode goes exponentially to zero at 
X —?> ±00 




exp[ua( 5 )x] , X ^-00 
exp [—Uf,((?)a:] , x —>■ -l-oo 


(47) 


As in the bosonic case, for what concerns the bound states relative to the vacuum | a), the only poles 
that matter are those coming from the behavior of this function (relative to the kink of lowest mass) at 
X —)■ — 00 : in particular, the first pole of G^f^{9) is simply fixed by the asymptotic value Va 


9i 


in 



and it is in the physical strip if 


Vaig) 


Vaig) 


< 1 . 


When this happens, we have two Q = 0 fermionic bound states of mass 


mf{g) 


± 2 M:, sin 



(48) 


(49) 


(50) 


Few comments on this formula: 


1. The ± signs refer to the paired solutions of energy ±E. From now on for simplicity we will concen¬ 
trate the attention only on one of them, say the positive energy bound state. 


2. When rja < 1, we have that | m{ (g) \< 2M*j, i.e. we are in presence of genuine Q = 0 fermionic 
bound states made of the kink-antikink of lowest mass carrying Majorana fermion number. 


3. In the very deep semi-classical limit, when —)■ 00 , expanding in series the formula (50), we have 

that I mi \— Va, which is of course the mass usually associated to the elementary fermion '0 in the 
(5 = 0 sector at the vacuum | a). 

4. If, on the contrary, moving the coupling constant g, we reach a value gc for which rja^gc) > 1) then 

I mi I overpasses the threshold value and leaves the physical spectrum. Therefore, similarly 

the purely bosonic case [2H], the picture is as follows: the particle associate to the fermion present 
in the Lagrangian (26) has to be considered, on a general ground, not an elementary particle but a 


bound state of the two kinks with different fermion parity, and it belongs to the physical spectrum 
of the theory if and only if the condition (49) is satisfied. 
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b> 

a) 


r~L 


KtbK^a) 


b> 

a) 




b> 

a-r) 

. _r" 

r 



Kba) 

a_) 

b> 


“C. 

a-r) 




FIG. 9: Two-particle fermionic kinks in the graphical representation of the original vacua (figures on the left) and 
of the unfolded vacua (figures on the right). 


In order to get the other poles eventually present in the Fourier transform of (x), we need to analyse 
its exponential approach at x —>■ —oo. This behavior evidently depends on the interaction term V(<p(x)) 
and, in absence of any further information on this function, we cannot reach any general conclusion. 
However, we can pin down the general pattern of the asymptotic behavior of the zero-mode in the very 
important cases when V((p) admits a series expansion around the vacuum values Indeed, assume 
that nearby the vacuum values the interaction term can be expanded as 


VabiTi^)) = Va 


■ Kbivix) - Ta'’) + TiVabiTix) - Ta'’f + 


(51) 


Substituting now in this expression the quantity {v?(x) — with the exponential approach given in ([t), 
we see that all terms but the first are expressed in terms of purely exponential terms which are multip e 
of a single one. Namely, the expression of Vab(<p(x)) around x -T —oo can be generally written as 


Vab{T{x)) = -Cg -f , X-oo , (52) 

n—1 


are coefficients determined 

(a) 


where uia is the curvature of the bosonic potential U{ip) at ip = ip^\ while d„ 
by combining both the various derivatives of the interaction term and the expansion coefficients pn 
in Q. However, we need an extra step to control the fermion zero mode itself, namely we need to evaluate 
the integral of Vab(.T) for x -T —oo. This can be easily done using eq. (52): indeed, disregarding the various 
constants of integration and introducing other appropriate coefficients the general expression of this 
integral can be cast as 


J Vab{t)dt 

Xq 


-VaX + Y, 
n=l 


X —>• —CX) , 


(53) 


where all terms, but the first one, exponentially small. Therefore, the final result is as follows: when 
Vabi^) is analytic around the vacuum configurations ip^\ we have the following asymptotic expansion 
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for X —)■ —00 of the fermionic zero mode 


exp 


= Atp e'’ 


~ X 


- 

1 

1 

1 

= A ij} exp 

-VaX + Y, 
n—1 


1 + t„ 


(54) 

(55) 


n—1 


where are the final coefficients resulting from the series expansion of the exponential function and the 
coefficients dn- The presence of these exponential terms induces poles in the function localised at 


= iTr{l -Tja-n^a) ■ 


(56) 


They are on the physical strip as far as n < (1 — ?7a)/faj whose integer part is therefore the total number 
of fermionic bound state in the Q = 0 sector 


= 


1 - 




(57) 


The general semi-classical expression of the masses of the fermionic bound states is thus given by 

n = , (58) 


which can be also written as 


< = 2*4Vi. . 

2 


= rui 


sin(d^) 


n = l,2,...,N^ 


(59) 


This is the main result of the paper. Comparing this mass formula with the one of the bosonic bound 
states, eq. (211, we see that the only significant difference is the shift induced by the quantity rja- In the 
next sections we are going to discuss a series of interesting QFT and study the various patterns of bound 
states that emerge in those cases. 


VI. SYMMETRIC WELLS 


The first theory we want to discuss is the familiar symmetric double well bosonic potential that conve¬ 
niently we choose to write as 


U{ip) 



(60) 


We denote with | ±1) the vacua corresponding to the classical minima = ±i^ = The 

kink/anti-kink solutions interpolating between them are 


P-a,a{x) = a (fi tanh 

where we have chosen the integration constant Xg = 
logical configurations is given by 


mx 

.71. 


a = ±1 


(61) 


0. Using eq. (12), the classical mass of these topo- 


Mo 



(62) 


Keeping into account the finite quantum correction of this quantity [251 [3l], the kink mass becomes 


M 


•\/2 

“ 3 “ 





(63) 



















Conveniently defining 
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c = 


27r 4^3 


and the dimensionless quantities 


a = 


3A2 


Trm^ 


the mass of the kink can be finally expressed as 

M = 


C = 


\/ 2 m 

TT^ 


>0 , 


1 — TTCg ’ 


(64) 


(65) 


In this example the kink and the anti-kink solutions are equal functions (up to a sign), and therefore 
their Fourier transforms have the same poles. Hence, the spectrum of the neutral particles will be the 
same on both vacua | ±1), in agreement with the Z 2 symmetry of the model, and is given by 


TO® „ = 2M sin 


TT^ 


n = 1, 2,... 7V|‘ 


( 66 ) 


The interested reader can find a more detailed analysis of this bosonic sector in [28j . 

Let’s now add the fermion and compute its spectrum (for a previous semi-classical study of this problem, 
close to the spirit of the DHN approach [5^, see [22] )■ Since V{ip) - evaluated on the kink solution - must 
assume opposite values at a; —)■ ± 00 , the simplest way to satisfy this requirement consists of a potential 
simply proportional to the bosonic field ^p{x) itself 

V{ip) = g(p{x) , (67) 

where g is another coupling constant. So, for the asymptotic values of H at a: —)■ ±00 we have 

Va = -V-a = gv ■ ( 68 ) 

Let’s now briefly analysed the fermionic spectrum in the Q — 0 and Q = 1 sectors. 

Fermionic bound states in Q = 0 sector. Given V{(p), we can compute the fermionic zero mode, 
which in this case can be expressed in a closed simple form 


V’-i.a(a’) = M 


1 


cosh^ 


where r is the ratio of the couplings 


_ a 
A 


Posed ijJa = V2m, its behaviour at a: —)■ —00 is given by 

~ 2’’ H [1 - r -k r(r -k 1) -k • • • ] 

Hence the poles of the matrix element G-aa{9) are located at 

6»„ = i7r(l-?7-n^) , 

with 

^ ttM ’ ttM ’ 

and the masses of the fermionic bound states at both vacua | ±) are given by 

'TT{r] -k n^)' 


TO^ „ = 2M sin 


n = l,2,...,iV^ 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 
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Comparing this spectrum with the one of bosons, eq. ( 661 , we see that the masses of bosons and fermions, 
as well as their number, are in general different, see Figure |10[ There is however a special case where 
the two spectra match exactly: it is when rj = ^, alias when g = X. As we will see in Section |VIII[ this 
is not a coincidence but it turns out to be the condition under which the theory is invariant under a 
supersymmetric transformation. 

Fermionic bound states in Q = ±1 sector. In the case of the symmetric wells it is possible to solve in 
a close form the set of Schrodinger equations (42) and hnd the fermion spectrum in the first non-trivial 
Q = ±1 topological sectors. The computation can be found in the original paper |25) or in a more recent 
one [35j : in the latter, in particular, to get the spectrum of the Dirac Hamiltonian it has been used 
the technique of shape-invariant potentials. Here we simply report the final results for the bound state 
energies which depends upon only on the ratio r of the two coupling constants 


fc = 1,2,... [r] 


(75) 


Ek = \/2rk — 

v2 

The upper value of these energies coincides with Emax 
so the entire discrete spectrum is in the interval {—Emax, Emax) as shown in Figure 11 Notice that such 
bound states only exist when r > 1, otherwise they are absent. In particular, there is none of them at 
the supersymmetric point g = X. 


= Vai which is the effective mass of the fermion, 


VII. ASYMMETRIC WELLS 

It is interesting to see what happens when the bosonic potential has two asymmetric wells | a) and | b): 
in this case the asymptotic behaviors of the kink at a: ± oo are different and therefore one should expect 
to find two different spectra piling up on the two vacua. For the purely bosonic model this case has been 
discussed in [28], where it has been shown in particular that there could exist a range of the coupling 
constant where one vacuum (say | a) ) has no bound states on top of it, while the other vacuum | b) 
can have instead one or more. Such a situation is realised, for instance, in the Tricritical Ising Model, 
once perturbed away from criticality by the sub-leading magnetic field [861137j . Here we analyse such an 
asymmetric potential when there are also fermions. Let’s first briefly remind the phenomenology of the 
purely bosonic case. 

As bosonic potential U{(p) of this model we choose 

U{^) = y (^ + «x)' (t’ - + Cy) . (76) 


E 


(a) 


(b) 


FIG. 10: Spectrum of boson and fermion bound states in the Q = 0 sector for: (a) generic values of the couplings 
g and X; (b) for g — X 
























FIG. 11: Discrete spectrum of fermions in the Q — 1 sector, including the zero mode for g/X = 6 


Rescaling the field as ^p{x) —>■ we can bring the Lagrangian in the form 


£ = 


A4 




(77) 


The minima of U{(p) are localised at = —a and pf'’ = b and the corresponding ground states will 
be denoted by | a) and | b). The curvature of the potential at these points is given by 


U"{—a) = Wq = (a + &)^(a^ + c) ; 
U"{b) = ujl = {a + hf{b‘^+c) , 


(78) 


and therefore for a ^ b, we have two asymmetric wells, as shown in Figure [I^ From now on we choose 
the curvature at the vacuum | a) to be higher that the one at the vacuum | b), i.e. a condition realised 
when a > b. 

The kink equation is given in this case by 


^ = ±{p + a){p - b) + c . 


(79) 


The kink interpolates between the values —a (at x = —oo) and b (at x = +oo) while the anti-kink does 
the viceversa. However, the behaviour of the anti-kink at x = —oo is different from the one of the kink, 
since the two vacua are approached differently: indeed, Pab(x) approaches the vacuum | a) as 


Pabix) ~ -a -I- riexp(a;aa:) 

while the anti-kink approaches instead the vacuum | b) as 

Pbaix) ^b- Bexp{LObx) 


X —>■ —oo 


X —>■ —oo 


(80) 


(81) 



FIG. 12: Example of p^ potential with two asymmetric wells and a bound state only on one of them. 
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where A and B are two positive constants. Since wq 7 ^ wi, the asymptotic behaviour of the two solutions 
ends up in the following poles in their Fourier transform 


■^(<Pab) 


■B(ipba) 


ik 


(82) 


-B 


uJh + ik 

A dimensional analysis fixes the mass of the kink to be 


"= v" 


(83) 


where a (that depends upon a^b,c) can be computed taking the integral (12). Using this expression of 
the mass and expressing the poles given above in the variable 9 we have 


qA) ~ 


m \ 


ttM / 


= ITT 1 — UJa 




arrr 




(84) 


Therefore, if we choose the coupling constant in the range 

TTCX 


< 


< 


7ra 




(85) 


the first pole will be out of the physical sheet whereas the second will still remain inside it. Namely, the 
theory will have only one neutral bosonic bound state, localised at the vacuum | b). Said differently, an 
antikink-kink configuration produces a bound state whereas a kink-antikink does not. 

Let’s now consider how the fermion may change the previous picture. As potential term V ((/?) we 
consider once again for simplicity U(i^) = gip{x), so that 


Va = -V{-a) = ga’f , 
V, = V{b) = gb^ . 


( 86 ) 


This time there will be a zero-mode corresponding to the kink configuration ipab(x) and another one 

" 0 ^°^ relative to the antikink (pba(x), with the behavior at ac —>■ =Foo of one equal to that one of the other 
at a; —)■ ± 00 , so that 


'^ab(^) = Alp exp 


'>Pba(^) = ^'0 exp 


- J Vabit)dt 

L Xo 

X 

- J Vba{t)dt 

L Xo 




Alp e" 


1 + tn ' 
V n=l 

1 + tn t 


— 00 


(87) 


( 88 ) 


Therefore, posing gk = Vk/i'xM) and = ^kK.'xM) {k = a,b), we find two different fermion spectra on 
the two vacua 


(rik) = 2M sin + nk^k)^ , rife = 1 , 2 ,... 


1 - r?fc 
■Cfc 


(89) 


Varying the coupling constants A and g one can get various and also curious situations at the two vacua: 
for instance, imagine that we are in the condition (85), where there is no bosonic bound state on the 
vacuum | a) and only one on the vacuum | b) . But we can easily revert the situation for what concerns the 
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FIG. 13: Asymmetric potential, with one fermion bound state at the vacuum \ a) and no boson bound state, and 
a boson bound state at the vacuum \ b) but no fermion one. 


fermionic spectrum! Alias, if we choose the ratio of the coupling constants g and A such that it satisfies 
the condition (where C = Tram'^(X^) 


C 

b 


9 C 

< T < - 
A a 


(90) 


then the vacuum | a) will have a fermionic bound state while the vacuum | b) has none, as shown in 
Figure [T^ 


VIII. SUPERSYMMETRIC THEORY 


Let’s consider now QFT’s which are invariant under a IV = 1 supersymmetry transformation. In this case, 
the bosonic and fermionic fields, together with a real auxiliar field F(x) can be conveniently organised 
into a real superfield $(x, 0 ) that admits the expansion 


$(a;,0) = ip{x) + 9tjj{x) + ^0 9F{x) 


(91) 


The space coordinates x^ = (a:°,a;^) and the two real Grassmann coordinates 0a = {0i,02) describe the 
iV = 1 superspace. Under a translation in superspace 


X^ ^ X^ + i^J^0 , 0a ^ t 

the variation of the superfield is given by 

5$(a::,6») = ^a Qa^ix,0) 


a F ^a 


(92) 


(93) 


with Qa = dld0a + ii'j'^0)a du- most general action invariant under the supersymmetric transfor¬ 
mation (93) is given by 


A = / (Fx(F0 


{Da^)Da<^ + W(<^) 


(94) 


where f (P900 = 2, with the covariant derivative Da given by 




(95) 


bF($) is the so-called superpotential, that we assume to be an analytic function of $. Integrating on the 
Grassmann variables, one arrives to the following expression of the action 


A = y ^^ 2 :1 i F FF"^] +FW {(p) - ^W” 


(96) 
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where W ((p) = dW((p)/dip, etc. Finally, elimitating the auxiliary field F from its algebraic equation of 
motion, i.e. substituting F —> —W ((p) in the above expression, and rescaling for convenience the fermion 
field as 'i/' —> it yields the general form of the lagrangian density for a supersymmetric theory given 

by 




Associated to the transformation (92) there is the conserved supercurrent 


J^{x) = 


and the conserved supercharges 


1 ro 


Qa = dx J, 


(97) 


(98) 


(99) 


Together with the stress-energy tensor 

r^^(x) = 1[{do^ipf - ^2] , (100) 

and the topological current 

= -e^^^Fd^ip = e^^''d^W{ip) . ( 101 ) 

they close the supersymmetry algebra [38j 

{Qa, 0/3} = 2(7A)a/3 + 2z(75)a/3 ^ , (102) 

where = f x) dx^ is the total energy and momentum, whereas Z plays now the role of the 
topological charge 


2 ab = Ie{^)dx^ = [vb( 7 ^)]^:: = w{p,k)-w{ipa ). 


(103) 


Although closely related to the topological charge introduced in (151, depending on the model Z may 
however differs from it. We will comment on this fact on the models considered later. A convenient and 
explicit form way of the supersymmetry algebra (102) is the following 

Ql = P+ , Ql = P- , Q+ Q- + Q_ Q+ = 2Z , 
where we have used the notation (Qi, Q 2 ) = (Q_, Q+) and P± = Pq ± Pi- 


(104) 


IX. GENERAL RESULTS IN SUSY THEORIES 


The existence of an algebra that relates the bosonic and fermionic sectors poses of course important 
constraints on the theory. Here we briefly collect some of these consequences. First of all, comparing the 
Lagrangian eq. (971 with the original one (p^, we see that in A = 1 SUSY model both potentials U{ip) 
and V{<p) come from the same function W((p), alias 


U(p.) = , 

V(<p) = W"(<p) . 


(105) 


U((p) is then intrinsicaly positive and it is well known its relation with the spontaneously breaking of 
supersymmetry |33) : in fact, when I7((p) has zeros, they are the true vacua of the theory and supersym¬ 
metry is unbroken; viceversa, if U((p) has local minima that are not zeros of this function, supersymmetry 
will be spontaneously broken at these minima. The local minima may be regarded as meta-stable vacua, 
as far as it exists a true vacuum somewhere in the landscape of U{ip). 

Focusing now the attention on the bosonic and fermionic bound states in Q = 0 sector, let’s show that 
the general pattern of their spectrum is fixed by the following two results. 
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1. An identity concerning the potentials. For all SUSY models it is easy to prove that for the 
fermionic and bosonic potentials we have 

Va = ±UJa ■ (106) 

Such an identity is easy to prove since 

Va = Via) = W"ia) , (107) 


while 


Wa = 0 = iW"ia))^ + W'ia)W'"ia) = (lU"(a))2 , (108) 

where the last step comes from W'(a) = 0, since (pa is one of the vacua of W’ip). 

2. Exact value of the classical mass of kinks. In SUSY theories the classical mass of the kinks can 
be exactly computed entirely in terms of the algebra [38] . Indeed the minima of U{(p) are connected 
by supersymmetric kinks, solutions of the so-called BPS equations 


(Q+±Q_) \ Kab) = 0 , 


(109) 


where the ± refers to the kink | Kab) and the antikink | Kba) respectively. The bosonic component 
of these equations is just the familiar condition (§ 


dpab 

dx 


iF{(pab) = ±W ipab) ■ 


Thanks to the identity 


P+ + P_ = iQ+±Q_)^T2Zab 


( 110 ) 


( 111 ) 


that follows from the supersymmetry algebra (104), the classical mass 
simply expressed by their topological charge: for the kink (anitkink) at 
given by P+ + P- = 2M, whereas the right hand side, using eq. (|109[), is 


Mab of the kinks is then 
rest, the left hand side is 
equal to ^2Zab, therefore 


Mab = \Z, 


ab \ 


( 112 ) 


Notice that this expression is a particular case of the general formula (12), given that i n sup ersym- 
metry yj2Uip) is equal to Wip) and the topological charge Zab is expressed by eq. (103). Let’s 
mention that it has been long disputed whether this result - true at the classical level - remains 
still valid once included quantum corrections (the interested reader may find a vivid summary of 
the story of this problem and the relevant references in the Chapter 3 of the book [35]). The final 
word on this problem seems to be as follows: at the quantum level, the mass of the kinks gets 
corrected Mab Mab + SMab but nevertheless satishes an equation as (112)) if one takes into a 
proper account that the superpotential W (and then Z) gets a correction too. 


Degeneracy of the bosonic and fermionic bound states. We can now used the two results given 
above to reach a very general conclusion for the spectra of supersymmetry theories. Namely, given the 
identity Wa =| fa | and the expression (112) of the mass of the kinks, for SUSY theories we have 


Va — Ca 


(113) 


and therefore the spectra (21) and (57) of bosonic and fermionic bound states at each true ground state 


of zero energy are always necessarily paired. 
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FIG. 14: Structure and adjacency diagram of the vacua in the SUSY Sine-Gordon in the periodic interval of the 
model. This structure has to be periodically repeated for the other vacua. 


X. INTEGRABLE SUSY MODELS 


There is a set ot theories where we can test our results: those are the integrable exactly solvable model. 
Here we focus the attention on a paradigmatic example of these theories, alias the SUSY Sine-Gordon 
model, whose superpotential W ($) 


cos(A$) 


gives rise to the SUSY Sine-Gordon Lagrangian 

n2 


-C = ^ ^ sin^ X(p -|- - rmpif cos X(p 


(114) 


(115) 


A big deal is known about this model, in many different contexts (see, for instance, [ISHU): in particular, 
the issue of the finite correction to the kink mass ( |112[) h as been clarified in [46] while the exact S'-matrix 
and the related spectrum have been determined in m So, the exact spectrum of the Q = 0 sector of 
this model is given by 


m^{n) = m^{n) = 2M sin 


(nXA 

, n = 1, 2,... 

■ 1 ■ 

-TT- 

— 

V 2 i 


.AA 


(116) 


where M is the mass of the kink while A^, the so-called renormalized coupling constant, is given by 


A2 


A^ = — 


1 




(117) 


Let’s see how the spectrum (116) is recovered by our semi-classical approach and, in the way to show 
this, also learning about other properties of this model which will be useful in the next section when we 
study the breaking of its integrability. 

The minima of the scalar potential U{g)) are localised at = kir/X. Si nce they are also the zeros 
of U{ip), they are supersymmetric vacua of the quantum theory. The theory (115) is invariant under the 
Z 2 parity p ^ —p and under the shift 


p ^ p + 


27rn 


(118) 


so we restrict our attention to the interval (0, 27r/A) of the field p. The elementary kink and the anti-kink 
solutions of ([^ of the SSG model are explicitly given by 


pf{x) = - arctan(e±’”’’’) 

A 


(119) 
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and connect the adjacent vacua = 0 and = tt/A, all other kinks or antikinks of the model being 
equivalent to them. In this model Z differs from the usual topological charge Q, since the latter counts 
the kink number while the former counts the kink number only modulo two |38j . 

In addition to the explicit expression of the kinks, there is also a close formula also the fermionic zero 
mode, since the integral of V{(pci{x) can be dealt as 


J V{^^'^{x)dx = J ^(V^) d.x = J = log [sin((^'='(a;))] 


( 120 ) 


V>‘^‘{xo) 


ip‘=‘{xo) 


where we have used eq.(llO) for the derivative of the kink solution and choose Xq such that ipci(xo) = 
7r/(2A). In this way, the zero mode can be exactly expressed as 


= Ai/jq sin [2arctan(e™'^)] = A ipo 


1 


cosh(TOa;) 


( 121 ) 


and, as a matter of fact, is equal to the zero-mode of the double well potential at the supersymmetric 
point g = X. 

About the nature of the SUSY vacua, we can take every vacuum = 2n tt/A to be non-degenerate 
(these are the vacua where V((p 2 n} > 0)j while we can unfold those located at = (2n -|- 1) tt/A 

(where < 0): the adjacency properties of the new vacua are shown in Figure 14 This way of 

unfolding the vacua is similar to what has been done in the Supersymmetry deformation of the Tricritical 
Ising Model |35j. The classical mass of these kinks, according to the formula (112), is given by 


M,i - . 


In references 


it has been computed the finite quantum correction of this mass 
Mci —i-M = 2m 

an expression that can be written as 


1 

A2 


1 

47r 


M = 


ttA^ 


( 122 ) 


(123) 


(124) 


where A is the effective coupling constant of the SUSY Sine-Gordon model given in (117). This formula 
closes now the circle, since all the minima of U(ip) has a curvature equal to uJa = m and therefore for the 
semiclassical parameter r/a and we have 


m 


“ ttM " 


(125) 


Hence, for this exac tly solvable model, the semi-classical spectra (21) and (74) precisely coincide with 
the exact ones (116). 


XI. NON-INTEGRABLE SUSY MODELS 

The SUSY Sine-Gordon model has provided an important check of our formulas for the mass spectrum 
of bosons and fermions. However the actual usefulness of the semi-classical formulas is of course in the 
absence of an exact solution of the model. This is, for instance, the case of the double well theory 
analysed in Section |vg even though when g = X this model can be elegantly put in a SUSY form, with 
the super-potential 

9 

1F(4>) = A$3 - —4> , (126) 

2 a 

it remains nevertheless non-integrable and therefore the discussion done in Section [VI| becomes particu¬ 
larly valuable. 
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(a) 


(b) 


FIG. 15: V{ip,p) for an irrational value of uj (here the golden ratio) and different values of the deformation 
parameter p: (a) p — 0.3; (b) p = 0.7. 


While the model defined by eq. (1261 has no parameter and it is simply non-integrable, here we are 
interested to study a class of non-integrable SUSY models defined by a Lagrangian C{p) depending on a 
parameter p that varies on the interval (0,1), such that the corresponding theory is integrable at both 
extrema, p = 0 and p = 1, but non-integrable otherwise. A representative of this type of models is the 
multi-frequency SUSY Sine-Gordon model m. with super-potential given by 


W{^,p) = m [ —^ cos(A<i)) -I- ^ cos{P^) 


(127) 


The associated bosonic potential is given by 

/1-p 


U(^,P) = ^ 


cos(A$) + ^ sin(^$) 
A p 


(128) 


We have chosen a proper parameterization in such a way that the curvature at the origin is always given by 
m^. In the flow of p from 0 to 1 one expects that the spectra at the different vacua move correspondingly 
and we are interested to see whether we can reach a certain control of this situation. 

First of all it is obvious that a relevant parameter is the ratio of the two frequencies 


P 

"= A 


(129) 


since we can a lways rescale the bosonic field as i p — » (p/X and manage to have lu as frequency of the second 
term in (1281. If w is irrational, the potential ( 128| ) loses any periodicity in ip and, apart the residue Z 2 
summetry (p —^ — (p, it generically gives rise to an infinite number of vacua that may move in a complicate 
way by varying the parameter p (see, for instance. Figure 15). Such a situation is the more complicated 
one but, using the continue fraction representation of any real number in terms of rational numbers, it 
can be studied as a limit case of the rational ratio of the frequencies, a situation to which we turn our 
attention. 

Let’s consider then the case when w is a rational number, uj = ^ (say with q > p). In this case, the 
periodicity of the original potential stretches to 


I = 



(130) 


As noticed in m, an important characteristic of the field theory described by ( |128[ ) is that the kinks 
initially present at p = 0 also persist for finite values of p, regardless of the ratio ui of the two frequencies 
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- a result that is strikingly different of what happens in the purely bosonic multi-frequency Sine-Gordon 
model [18j . The proof of this feature of the model, based on its Supersymmetry and pursued through the 
Form-Factor Perturbation Theory, can be found in |16] . An intuitive way to understand this robustness 
of the kinks is in terms of the persistence of the original zeros of the bosonic potential U{if,p — 0) for 
finite values of p. Moreover, as discussed below, when they disappear they do so in pairs and this happens 
only at some critical values pm of the parameter p. This observation helps in understanding the evolution 
of the kinks and the number of stable particles at any given p. To make it quantitative, call A4q{p) the 
number of zeros of U{ip, p) at a given value of p in the extended interval of periodicity (130). Clearly 


Mo{Q)=2q+l 


Afo(l) =2p+l , 


(131) 


so that from the beginning to the end of the flow there has been a total variation AA4q = 2 {q — p) of the 
zeros. Since the existence of the kinks is linked to the zeros of U{ip,p), a variation of the latter implies 
the disappearing of some of the former along the flow from p = 0 to p = 1. Being the kinks topological 
excitations, their disappearance will signal the occurrence of a phase transition at some critical values of 
the parameter, p = pi™\ which are the values where the number of zeros jumps by a step of 2 units. In 
fact, two zeros collide and then move on imaginary values; the barrier between these zeros vanishes and 
correspondingly the kink/antikink that connect the two colliding vacua become massless. 

The persistence of the kinks of the original theory at weak coupling does not mean that their mass 
remain untouched. In particular, the mass degeneracy of all the kinks of the integrable SUSY Sine-Gordon 
model at p = 0 is split when p is switched on, and we will have long-kinks (overpassing higher barriers) 
and short kinks (overpassing lower barriers), as shown for instance in Figure 9. The actual computation 
of the new classical masses of the kinks can be done by employing their topological charges, i.e. by using 
eqs. (103) and (112), together with the positions of the new zeros of U{(p,p). 

We can now write down the steps to take in order to see how the spectrum evolves by varying p: 


1. Determine in the interval (1301 the zeros p- 
are associated to the stab( 


(p) (n = 0,1,..., fc) of U(ip, p) by moving p. Such zeros 
e vacua | n) at that value of p. Their number depends on p and changes 


discontinuously of 2 units at some critical values pm . Notice that the zeros placed at = 0, 


.( 0 ) 


(fig ' = 7 rq/A and (p2q^ = ^'KqjA are always there, for any value of p. 

2. Compute the curvatures ujn{p) of the potential U{(p,p) at these zeros. 

3. Use the zeros previously determined to compute the classical masses Mat of the kinks connecting 
the various vacua according to the topological formula (112) 


Mab = \W{p'^^^)-W{pf^) 


(132) 


The quantum corrections may change these values but in the semi-classical approximation (A —> 0) 
one shall expect these corrections to be small. 

4. For each vacuum | a), identify the kink pab with the lowest mass M*;,. 

5. Use the two previous data to compute for each vacuum | a) the quantity fa = t^a/(7’‘M*). 


6. The number of stable bosonic and fermionic particles at each vacuum | a) is then given by Na = 


and their mass is obtained by using the semi-classical formula (21) or (74). 


7. Identify the critical values pm when two zeros collide and disappear afterward. At these values 
the related kinks become massless and this phenomenon produces a singularity in certain ^a’s, as 
discussed in more detail in the next Section. 

It is interesting to see how this protocol determine the evolution of the spectrum in two simple but 
significant examples. 


XII. VACUA REALIZATION OF SUSY 

The disappearence of the vacua moving the parameter p poses the question how SUSY is realised on the 
various vacua and whether it can get spontaneously broken in some of them even though it was originally 
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FIG. 16: Evolution of the potential p) for the SUSY Double Sine-Gordon model: (a) p = 0.3; (b) p = 0.5 and 
p = 0.8. In the figures there are drawn for simplicity only the kinks, in red the long-kinks, in blue the short-kinks 


exact on all the initial vacua. The answer is generically afFermative and this because, once two zeros 
collide at the position at some pc, for p = pc + e we have V , pc + e) > 0. This means that SUSY 
will be spontaneously broken at this point of the field space and, in the presence of other zeros, the new 
minimum at becomes a meta-stable vacuum state. There is a notable exception to this scenario and 
it is related to the arithmetic properties of the ratio lo |16j . Namely, when (g — p) is an odd number, 
among the pairs of colliding zeros there will always be a couple of them placed just on the right and on 
the left of Lp'f'^ = nq/X: when these zeros collide at some critical value Pc they strangle the zero at 7rg/A 
that is in between, so that ipq^'^ remains a zero even after their collision. This situation does not happen 
instead when (q — p) is an even number. 

Since in the ultraviolet, all multiple SUSY Sine-Gordon models may be regarded as perturbation of 
the Superconformal Field Theory with central charge c = 3/2, it was argued in [16]) that these theories 
give rise in general to a sequence of phase transitions that locally recalls the spontaneously symmetry 
breaking that occurs in the Tricritical Ising Model: moreover, when (g — p) is an odd number, there is 
also an extra phase transition similar to the one of the gaussian model. 

Let’s now study in a certain detail two models which are somehow representative of the two classes 
where (g — p) is an odd or an even number. 


A. Double Sine-Gordon Model 

The first example is the Double SUSY Sine-Gordon model, with a bosonic potential given by 

^ ^^-y^sinA(p-k y sin^^ . (133) 

In this case w = ^ and {q — p) = 1. As we will see soon, this model has only one critical value of the 
parameter p at pc = 1/2. Three different snapshot of the potential taken at p < pc, p = Pc and p > pc, 
shown in Figure [Tbl help in understanding and visualising the evolution of the model: 
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FIG. 17: Root relative to the vacuum \ 1) of the SUSY Double Sine-Gordon as a function of p. For p > pc = \; 
this root coincides with . 


a When p < pc, in the enlarged interval of periodicity (130) there are 4 stable vacua, identihed as in 


Figure 16 a. The vacua | 0) and | 1) are connected by long-kink/antikinks | and | The 


vacua I T) and | 2), as well as | 2) and | 3), are connected instead by the short-kink/antikinks | 

I and | K 23 ), \ K^^)- In this model, the reflection symmetry with respect to the vacuum | 2 ) 

provides an obvious identihcation of the latter four kink s. Th e vacua | 0) and | 2) correspond to the 
two fixed zeros (/jq = 0 and ip^'^ = tt/A of the potential (133) for all values of p. 


b At p = Pc = 5 , the two vacua | 1) and | 3) collide, strangling between them the vacuum | 2 ). 
Correspondingly the kinks | and | AT^) (together with their anti-kinks) become simultaneously 
massless. The critical value pc is identified as the value for which the curvature at the vacuum 
I 2) vanishes. The vacuum | 2), according to our general analysis, survives the collision and the 
disappearance of the vacua | 1) and | 3). 

c When p > pc, there are only | 0) and | 2) as stable vacua in the theory, connected by a single kink 
I Kq2) and the corresponding antikink | K^q). 

Let’s see how the spectrum at each vacuum evolves by varying p, choosing a particular value of A, say 
A = 2 -^ 7 r/ll. It is worth saying that the choice of any other value of A does not change the overall 
picture but just its details, as will become soon clear from the discussion below. 

We need first to determine the position of the vacuum | 1) (and by symmetry with respect to 
the vacuu m | 2 ), the other zero This can be done numerically, finding the non-trivial zero of the 

potential (133): the result is plotted in Figure [l7| Once is known, using eq. (112) we can compute the 
masses of the long and short kinks and their behavior versus p is shown in Figure 18 Computing 


as well as the curvature at each vacua and always using the shortest kink for each vacuum, we can then 


determine the important parameters ^a(p) and the numbers Na{p) = 
vacua. Let’s now discuss what happens at each vacuum separately: 


of the bound states at each 


1. the vacuum | 0) exists for all values of p, therefore it makes sense to consider Nq{p) along the 
entire flow. As shown in Figure [T^ in this case No(p) is a increasing staircase function which jumps 



FIG. 18: Mass Ml (in red) and Ms (in blue) of the long and short kink respectively, as function of p. Notice that 
Ms = 0 for p> Pc. 
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FIG. 19: Number of bound states No and Ni as functions of p. 


discontinuously each time that p overpasses certain special values At the beginning of the flow, 
the number of bound states is fixed by the first frequency A and it is given by Nq{0) = 3, while at 
each a new bound state enters the spectrum coming from the continuum, so that at the end of 
the flow the total number of bound states is A'o(l) = 13, which is the number of bound state of the 
SUSY Sine-Gordon model relative to the second frequency 


2. the vacua | 1) and | 3) exist only for p < Pc and therefore it makes sense to talk about the number 
of the bound states A^i (p) on top of these vacua only for this range of p. For p > pc we may regard 
this number to be instead equal to zero. The result is the step-wise decreasing function also shown 
in Figure 19 which starts as before from A^i(O) = 3 and ends down then to zero. 


3. the number of bound states N 2 {p) on the vacua | 2) presents a discontinuity at p = Pc for the 
merging of three vacua at this value of p. For p < pc, this number N 2 {p) starts at p = 0 with the 
value 3 and decreases toward zero at pc, similarly to Ni{p). But at pc this number diverges, simply 
because: (a) the curvature of this vacuum at that value vanishes and (b) the mass of the kink that 
has to be employed to compute for p > pc is no longer Ms (that is zero) but Ml, which is instead 
finite. So, for p > pc, N 2 {p) decreases from values infinitely large toward the final value ^'2(1) = 13 
dictated by the second frequency of the model, see Figural^ 


With all these information one can then apply the semiclassical mass formulas to compute the actual 
masses of the bound states at each vacuum. Here we present the result only for the vacuum in the origin. 



0.0 0.2 0.4 0.6 0.8 1.0 


P 


FIG. 20: Number of bound states N 2 as function of p. Notice that N 2 diverges at p = pc. 

















FIG. 21: Evolution of the various masses of the bound states No at the vacuum in the origin versus p. In solid red 
line those which are always present in the flow, where in green dashed lines those corresponding to the additional 
bound states which start to appear only once p overpasses the special values p’^, even though in the Figure are 
plotted for all values of p. 


Figure 21 Notice that the three lowest lines are present all over the flow while there is a cascade of 


additional bound states coming from the continuum each time that p overpasses the values p* 


B. Triple SUSY Sine-Gordon Model 


This is the model (128) where the ratio of the frequency is w = 1/3. The landscape of the potential of 


this model for various values of p can be found in Figure One can see that, contrary to the previous 



10> |1> |2> |3> |4> 15> 

(a) 





csma 




|0> |1=2> |3> |4=S> 

(C) 


|0> |3> 

(d) 


a ma 



|0> |1> |2> |3> |4> |5> 

(b) 





FIG. 22: Evolution of the potential V{p,p) for the SUSY Triple Sine-Gordon model: (a) p = 0.02; (b) p — 0.05; 
(c) Pc = 0.2; (d)p = 0.4 and p = 0.8. In the figures there are drawn for simplicity only the kinks, in red the 
long-kinks, in blue the short-kinks 
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FIG. 23: Mass Ml of the long kink connecting the vacuum \ 0) to its nearest one. respectively, as function of p. 
Notice the discontinuity of the function at pc- 


model, here when two vacua collide at pc (| 1 ) and | 2 ), as well as | 4 ) and | 5 )), immediately after 
the potential is lift up, alias the two vacua disappear and nothing remains behind. This has the striking 
consequence that, not only the short kinks disappear at the critical point Pc but also that the long kink 
immediately after the transition is made of the two previous long kinks. This creates a discontinuity in 
the value of the mass of the lowest kink that has to be employed in the computation of No{p) and N 3 {p), 
as can be seen from FigureCorrespondingly, there is a jump of 2 units in the number of bound states 
that piled up on the vacua | 0 ) and | 3 ), see Figure 24 So, the breaking of SUSY at the colliding vacua 
I 1 ) and I 2 ) (and the companion | 4 ) and | 5 )) has an effect on the spectrum of the nearby vacua where 
SUSY is instead still exact. 


XIII. CONCLUSIONS 


In this paper we have derived the semi-classical formula (58) for the fermionic spectrum of an interacting 
QFT of boson and fermion. Together with the formula (21) for the bosons, these expressions provide 


interesting information on those quantum field theories made of degenerate vacua. Notice that the formula 
(58) remains valid also when the fermion is of Dirac type rather than Majorana, modifying correspondingly 


the expansion of the fermion in the presence of the kink, eq. (37), and using the more general relation 
between matrix elements given in eq. (|45|). So, these semi-classical formulas allow us to identify the particle 



FIG. 24: Number of bound states at \ 0) and \ 3). Notice the jump of 2 units at p^. 
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excitations and estimate one of their most important characteristics, i.e. their mass. Of course, in absence 
of integrability, not all the particles entering the semi-classical formulas are stable: to determine which 
of them will decay simply involves to identify the proper threshold of the decay channel - a task not 
particularly difficult but which we did not pursue here. The main difference with respect to the bosonic 
case analysed in [5H1 is the presence of the shift rja in the fermionic formula (58) with respect to the 
bosonic one (21), which may influence the value of the thresholds. However it is natural to think that, as 
in the bosonic case [55] , only few excitations will be generally stable on each vacua, all the others become 
resonances. 

With the fast advances of experimental works in quantum world - some of the recent achievements 
were briefly underlined in the Introduction - there may be soon the appealing perspective to realise 
many of the quantum field theories on a desk of atomic or material science laboratory. There has been 
in literature various proposals to address Majorana fermions in a context similar to the one discussed 
here, see for instance [30] for a proposal of observing Majorana bound states of Josephson vortices in 
topological superconductors, or m for the possibility to observe an emergent Supersymmetry from 
strongly interacting Majorana zero modes in the Tricritical Ising Model. If and when these or other 
proposals could be implemented experimentally , it would be extremely intriguing to check some of the 
features discussed here, among which: the patter of the bound states on the various vacua in a situation 
of local breaking of supersymmetry in models as the multiple Sine-Gordon theories; or, in a situation of 
asymmetric potential, the swapping of the bosonic and fermions excitations present on the two vacua. 


Acknowledgements 

I am very grateful to Michele Burrello for interesting discussions and for his careful reading of the 
manuscript. This work acknowledges the IRSES grant FP7-PEOPLE-2011-IRSES QICFT 295234. 




32 


[1] A.B. Zamolodchikov, Integrable field theory from conformal field theory, Adv. Stud. Pure Math. 19 (1989), 
641. 

[2] G. Delfino and G. Mussardo, The Spin spin correlation function in the two-dimensional Ising model in a 
magnetic field at T — Tc, Nucl.Phys. B455 (1995) 724-758 

(31 R. Coidea et al. Quantum Criticality in an Isinq Chain: Experimental Evidence for Emergent Ex Symmetry, 
Science, 327 (2010), 177-180 

[4] A.B. Zamolodchikov and ALB. Zamolodchikov, Factorized S matrices in two dimensions as the exact solutions 
of certain relativistic quantum field models, Ann. Phys. 120 (1979) 253. 

[5] M. Karowski and P. Weisz, Exact form factors in (l-hl)-dimensional field theoretic models with soliton be¬ 
haviour, Nucl.Phys. B 139 (1978) 455. 

[6] F.A. Smirnov, Form factors in completely integrable models of quantum field theory, Advanced series in 
mathematical physics, 14, World Scientific (1992). 

[7] T. Schweigler et al. On solving the quantum many-body problem, arXiv: 1505.03126 

[8] F. Halzen and A. Martin, Quarks & Leptons: An Introductory Course in Modern Particle Physics, John 
Wiley & Sons. (1984). 

[9] G. P. Engel, C. B. Lang, D. Mohler, A. Schafer, QCD with Two Light Dynamical Chirally Improved Quarks: 
Baryons, Phys. Rev. D 87 (2013) 7, 074504; 

S. Aoki et al.. Review of lattice results concerning low-energy particle physics, Eur.Phys.J. G74 (2014) 2890 
N. Brambilla et al., QCD and Strongly Coupled Cauge Theories: Challenges and Perspectives, Eur.Phys.J. 
G74 (2014) 10, 2981. 

[10] H. Bethe and E. Salpeter, A Relativistic Equation for Bound-State Problems, Phys. Rev. vol. 84 (1951), 1232. 

[11] G. Mussardo, Statistical Field Theory. An Introduction to Exactly Solved Models of Statistical Physics, Oxford 
University Press (2010). 

[12] G. Delfino, G. Mussardo and P. Simonetti, Nonintegrable quantum field theories as perturbations of certain 
integrable models, Nucl. Phys. B 473 (1996), 469. 

[13] G. Delfino and G. Mussardo, Nonintegrable aspects of the multifrequency Sine-Cordon model, Nucl. Phys. B 
516 (1998),675 

[14] D. Gontrozzi and G. Mussardo, Mass generation in perturbed massless integrable models, Phys. Lett. B 617 
(2005), 133. 

[15] G. Delfino, P. Grinza and G. Mussardo, Decay of particles above threshold in the Ising field theory with 
magnetic field, Nucl. Phys. B 737 (2006), 291. 

[16] G. Mussardo, Kink Confinement and Supersymmetry, JHEP 0708 (2007) 003. 

[17] V. P. Yurov and Al. B. Zamolodchikov, Truncated Conformal Space Approach to Scaling Lee-Yang Model, 
Int. J. Mod. Phys. A 05, 3221 (1990). 

[18] D. Gontrozzi and G. Mussardo, On the mass spectrum of the two-dimensional 0(3) sigma model with theta 
term, Phys. Rev. Lett. 92 (2004), 021601. 

[19] G. Delfino and P. Grinza, Confinement in the q-state Potts field theory, Nucl, Phys. B, vol. 791,(2007); 

L. Lepori, G. Toth, G. Delfino, Particle spectrum of the 3- state Potts field theory: A Numerical study, JSTAT 
vol. 0911:P11007. 

[20] P. Fonseca, A. Zamolodchikov, Ising spectroscopy. 1. Mesons at T < Tc, hep-th/0612304 
A. Zamolodchikov, Ising Spectroscopy II: Particles and poles at T > Tc, arXiv:1310.4821 

A. Zamolodchikov and 1. Ziyatdinov, Inelastic scattering and elastic amplitude in Ising field theory in a weak 
magnetic field at T > Tc: Perturbative analysis, Nucl.Phys. B849 (2011) 654-674. 

[21] S.B. Rutkevich, Baryon masses in the three-state Potts field theory in a weak magnetic field, JSTAT 1501 
(2015) 1, POlOlO; Two-kink bound states in the magnetically perturbed Potts field theory at T < Tc, J.Phys. 
43 (2010) 235004; Form Factor Perturbation Expansions and Confinement in the Ising Field Theory, J.Phys. 
A42 (2009) 304025. 

[22] G. Mussardo and G. Takacs, Effective potentials and kink spectra in non-integrable perturbed conformal field 
theories, J.Phys. A42 (2009) 304022 

[23] Z. Bajnok, L. Palla, G. Takacs, F. Wagner, Non-perturbative study of the two frequency sine-Gordon model, 
Nucl.Phys. B601 (2001) 503-538; 

B. Pozsgay and G. Takacs, Characterization of resonances using finite size effects Nucl.Phys. B748 (2006) 
485-523; 

M. Lencses, G. Takacs, Confinement in the q-state Potts model: an RG-TCSA study, arXiv: 1506.06477 

[24] G. ’t Hooft, Two-dimensional model for mesons, Nucl. Phys. B 75 (1974), 461-470; 

S. Huang, J.W. Negele and J. Polonyi, Lattice gauge theory and the structure of the vacuum and hadrons 
Nucl. Phys. B 307 (1988), 669-704; 

G. G. Gallan, N. Coote and D. J. Gross, Two-dimensional Yang-Mills theory: A model of quark confinement, 
Phys. Rev. D 13, 1649; 

W. Krauth and M. Staudacher, Nonintegrability of two-dimensional QCD, Phys. Lett. B 388 (1996), 808-812; 




33 


V. A. Fateev, S.L. Lukyanov and A. B. Zamolodchikov, On mass spectrum in ’t Hooft 2D model of mesons, 
J.Phys. A42 (2009) 304012. 

[25] R.F.Dashen, B.Hasslacher and A.Neveu, Nonperturbative methods and extended-hadron models in field theory. 
II. Two-dimensional models and extended hadrons, Phys. Rev. D 10 (1974) 4130; Particle spectrum in model 
field theories from semiclassical functional integral techniques, Phys. Rev. D 11 (1975) 3424; Semiclassical 
Bound States in an Asymptotically Free Theory, Phys.Rev. D12 (1975) 2443 

[26] J. Goldstone and R. Jackiw, Quantization of nonlinear waves, Phys.Rev. D 11 (1975) 1486. 

[27] G. Mussardo, V. Riva and G. Sotkov, Finite volume form-factors in semiclassical approximation, Nucl. Phys. 
B 670 (2003), 464. 

[28] G. Mussardo, Neutral bound states in kink-Uke theories, Nucl. Phys. B 779 (2007), 101-154. 

[29] R. Jackiw and C. Rebbi, SoUtons with fermion number 1/2, Phys. Rev. D 13 (1976), 3398. 

[30] E. Grosfeld and A. Stern, Observing Majorana bound states of Josephson vortices in topological superconduc¬ 
tors PNAS 108, 11810 (2011). 

[31] R.Rajaraman, SoUtons and instantons, Amsterdam, North Holland, 1982. 

[32] D. K. Campbell and Y. T. Liao, Semiclassical analysis of bound states in the two-dimensional a-model, Phys. 
Rev. D 14 (1976), 2093. 

[33] E. Witten, Dynamical breaking of supersymmetry, Nucl. Phys. B185 (1981) 513. 

[34] F. Cooper, A. Khare, U. Sukhatme, Supersymmetry and Quantum Mechanics, Phys. Rep. 251 (1995) 267-385. 

[35] F. Charmchi and S. S. Gousheh, Complete spectral analysis of the Jackiw-Rebbi model, including its zero 
mode, Phys. Rev. D 89, 025002 (2014). 

[36] M. Lassie, G. Mussardo and J.L. Gardy, The scaling region of the Tricritical Isinq Model in Two Dimensions, 
Nucl. Phys. B 348 (1991), 591. 

[37] F. Golomo, A. Koubek and G. Mussardo, The S-matrix of the Sub-leading Magnetic Deformation of the 
Tricritical Ising Model in Two Dimensions, Int. Journ. Mod. Phys. A 7 (1992), 5281. 

[38] E. Witten and D. Olive, Supersymmetry Algebras that include Topological Charges, Phys. Lett. B 76 B (1978), 
97. 

[39] M. Shifman and A. Yung, Supersymmetric SoUtons, Cambridge University Press, Cambridge 2009. 

[40] S. Sengupta, P. Majumdar, Conservation laws in the supersymmetric sine-Gordon model, Phys. Rev. D 33 
(1986), 3138. 

[41] C. Ahn, Complete S-Matrices of Supersymmetric Sine-Cordon Theory and Perturbed Superconformal Minimal 
Models, Nucl. Phys. B 354 (1991), 57; 

S. Shankar, E. Witten, S-matrix of supersymmetric cr model, Phys. Rev. D 17 (1978), 2134. 

[42] A.M. Tsvelik, Exact solution of (1-hi)-dimensional relativistic supersymmetric theories, Sov. J. Nucl. Phys. 
47 (1988), 172; 

A. Hegedus, F. Ravanini, J. Suzuki, Exact finite size spectrum in super sine- Gordon model, Nucl. Phys. B 
763 [FS] (2007), 330. 

[43] Z. Bajnok, C. Dunning, L. Palla, G. Takacs and F. Wagner, SUSY sine-Gordon theory as a perturbed conformal 
field theory and finite size effects, Nucl. Phys. B 679 [FS] (2004), 521. 

[44] M. Sakagami, Scattering of SoUtons with fermion number Nucl. Phys. B 207 (1982), 430. 

[45] A.B. Zamolodchikov, Fractional-spin integrals of motion in perturbed conformal field theory, in H. Guo, Z. 
Qiu, H.Tye (Eds), Fields, Strings and Quantum gravity, Gordon & Breach, 1984. 

[46] M.A. Shifman, A.I. Vainshtein, M.B. Voloshin, Anomaly and quantum corrections to solitons in two- 
dimensional theories with minimal supersymmetry, Phys. Rev. D 59 (1999), 045016; 

A. Losev, M. A. Shifman, and A. 1. Vainshtein, Counting supershort supermultiplets, Phys. Lett. B 522, 327 

( 2001 ) 

H. Nastase, M.A. Stephanov, P. van Nieuwenhuizen, A. Rebhan, Topological boundary conditions, the BPS 
bound, and elimination of ambiguities in the quantum mass of solitons, Nucl. Phys. B 5)2 (1999), 471; 

A. Goldhabera; A. Rebhanb, P. van Nieuwenhuizena, R. Wimmerc, Quantum corrections to mass and central 
charge of supersymmetric solitons, Phys. Reports 398, 179 (2004); 

A. Rebhan, P. van Nieuwenhuizen, R. Wimmer, The anomaly in the central charge of the supersymmetric 
kink from dimensional regularization and reduction, Nucl. Phys. B 648 (2003), 174; 

[47] A. Rahmani, X. Zhu, M. Franz, 1. Affleck, Emergent Supersymmetry from Strongly Interacting Majorana Zero 
Modes, arXiv:1504.05192 


